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We have proposed a novel numerical method to calculate accurately the physical quantities of 
the ground state with the tensor- network wave function in two dimensions. We determine the 
tensor network wavefunction by a projection approach which applies iteratively the Trotter-Suzuki 
decomposition of the projection operator and the singular value decomposition of matrix. The 
norm of the wavefunction and the expectation value of a physical observable are evaluated by a 
coarse grain renormalization group approach. Our method allows a tensor-network wavefunction 
with a high bond degree of freedom (such as D — 8) to be handled accurately and efficiently in the 
thermodynamic limit. For the Heisenberg model on a honeycomb lattice, our results for the ground 
state energy and the staggered magnetization agree well with those obtained by the quantum Monte 
Carlo and other approaches. 



The application of the density matrix renormaliza- 
tion group (DMRG) proposed by White [1] has achieved 
great success in one dimension in both zero and finite 
temperatures d, 0, 0|. However, in two dimensions, the 
application of the DMRG in both the real and momen- 
tum space [H, (sj has been limited only to small lattices. 
The error resulting from the DMRG truncation increases 
extremely fast with increasing size of lattice. To resolve 
this problem, the tensor- network state, which is an ex- 
tension of the matrix product in one dimension [7], was 
proposed[8j. In a tensor network state, a local tensor is 
directly entangled with the other local tensors on all the 
directions of the lattice. This leads to two problems in 
the treatment of the tensor network state. First, it is 
difficult to determine accurately all the elements of local 
tensors by any variational approach since the total degree 
of freedom of a local tensor increases exponentially with 
the dimension of the tensor. Second, it is even more dif- 
ficult to calculate the expectation value of any physical 
observable even if we know the expression of the tensor 
network wave function, since the number of summations 
over the basis configurations increases exponentially with 
lattice size. 

In this paper, we propose a novel method to handle 
the tensor-network wave function in two dimensions. We 
will show that the tensor network wavefunction of the 
ground state can be accurately determined by applying 
an iterative projection approach. This approach is sim- 
ilar to the time-evolving block decimation method that 
was used to determine the matrix product wave func- 
tion of the ground state in one dimension fioj Then we 
will generalize the classic coarse grain renormalization 
group approach proposed by Levin and Nave[ll| to the 
quantum system, and use it to calculate the norm of the 
wave function and the expectation value of any physical 
observable. This provides an accurate and efficient tool 
to determine the expectation values of physical quanti- 



tise from the tensor-network wavefunction of the ground 
state. 

Below, we will take the S=l/2 Heisenberg model on a 
honeycomb lattice 



h - Hj j , 



(1) 



h 13 = jSi • Sj - -h [(-ys^ + (-iys j:Z ] , (2) 

as an example to show how the method works. Here, (ij) 
means that i and j are the two nearest neighboring sites, 
h is the magnitude of a staggered magnetic field. It is 
straightforward to extend the method to other quantum 
lattice models with short range interactions in two di- 
mensions. We assume the tensor network state to have 
the following form 

I*) = Tr 11 X Xi X yi X Zi A 

XiUiZi [ m i]^XjyjZi [ m j] 

\mirrij). (3) 

A schematic representation of this tensor network state 
is shown in Fig. [TJ In Eq. ([3|), stands for the 

black/ white sublattice. vrti is the eigenvalue of Si Z . 
A Xi y iZi [mi] and B XjyjZj [mj] are the two three-indexed 
tensors defined on the black and white sublattices, re- 
spective. X a . (a = x, y, z) is a positive diagonal matrices 
(or vectors) of dimension D defined on the bond emitted 
from site i along the a direction. The subscripts a^, yi 
and Zi are the integer bond indices of dimension D (i.e. 
each running from 1 to D). A bond links two sites. The 
two bond indices denned from the two end points take 
the same values. For example, if the bond connecting i 
and j along the x direction, then X{ = Xj. The trace is 
to sum over all spin configurations {• • • , m^, m^, • • •} and 
over all bond indices. 
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FIG. 1: Schematic representation of a tensor network state 
on a honeycomb lattice. The lattice is divided into two sub- 
lattices, represented by the black and white dots, respectively. 
Each vertex, on which a spin state is inhibited, is connected 
with three neighboring vertices along three directions, labelled 
by x, y, and z. On each bond, there is a diagonal matrix (or a 
vector), A a , where the subscript a = x, y or z is a bond index 
of dimension D. At each vertex, a tensor representation of the 
spin state m, A x , y , z [m] for the black sublattice or B Xjy , z [m] 
for the white sublattice, is defined. A tensor network state is 
a product of all these bond vectors and vertex tensors. 



The ground state wavefunction can be determined by 
applying the projection operator exp(— tH) to an arbi- 
trary initial state In the limit (3 — > oo, exp(— rH)\^f) 
will converge to the ground state of H. However, this pro- 
jection cannot be done in a single step since the terms in 
H defined by Eq. ([2|) do not commute with each other. 
In real calculation, we will take a small r and apply this 
projection operator to \^f) iteratively for many times. 

Let us start by dividing the Hamiltonian into three 
parts 

H = H x + H y + H z , 

iGblack 

H a (a = x, y, z) contains all the interaction terms along 
the a-direction only. These terms commute with each 
other. From the Suzuki- Trotter formula, we can then 
express the projection operator as 



-Tff we -Tff, e -Tff Ve -Tir. +o(r 2 )- 



(4) 



This means that each iteration of projection can be done 
using exp(— rH a ) (a = x, y, z) in three separate steps. 

In the first step, the projection is done with H x . As 
only the two neighboring spins connected by horizontal 
bonds have interactions in H x , the resulting projected 
wavefunction can expressed as 



= Tr \\ J2 (rn / irn / j \e- Hi i T \m i m j ) 

i£b,j=i+x mirrij 

X Xi X Vi X Zi A XiVi z% [mi] B Xj y . Zj [rrij] \ m [m'j ) . (5) 
From this, a (D 2 d) x (D 2 d) matrix can be defined by 



J yiZi r m i ,yj Zjm\. 

(m • rrij \ e ~ Hij T \ <m % <m 3 ) 



mirrij x 



Xy i A-xyiZi 



i]X x B x . 



\KK, (6) 



where d = 2 is the total number of states of a S=l/2 
spin. Taking the singular value decomposition for this 
matrix, one can further express this S matrix as 



J yiZirrii,yjZjmj 



UyiZirniiX^xVxMjZjrrij (7) 



where U and V are two unitary matrices and A^ is a 
positive diagonal matrix of dimension D 2 d. 

Next we truncate the basis space by keeping only the 
D largest singular values of A^. Then we set the left \ x 
as the new A^ (x = 1 • • • D) and update the tensors A 
and B by the following formula 



A X y iZi [mi] = Xy^X^Uy, 

B xyjZj [mj\ = Xy^X-^V % 



yjZjr 



(8) 
(9) 



A flow chart of the above one-step renormalization of 
the wave function is shown in Fig. [2j The next two steps 
of projections can be similarly done with H y and H z , 
respectively. This completes one iteration of the projec- 
tion. By repeating this iteration proceduce many times, 
an accurate ground state wave function can then be pro- 
jected out. This iteration process is very efficient. The 
converging speed depends on the truncation error. In our 
calculation, we take r = 10 -3 initially and then gradu- 
ally reduce it to ~ 10 -5 to ensure the convergence of 
the wavefunction. The number of iterations used in our 
calculation is generally around 10 5 ~ 10 6 . 

Given the expectation value of a measurement 
quantity O is defined by 



(O) 



(10) 



We notice that both (Sfr\Sfr) and (*|0|#) are tensor- 
network functions. For example, 

/^|^r\=Tr TT T a , , ,T h , , , (11) 

\^|^/ -LJ- ^XiX\,yiy' v Ziz' i J -x 6 x , .,yjy' p Zjz' j -> \- LJ -J 

where the trace is to sum over all bond indices. Both T a 
and T b are D 2 x D 2 x D 2 tensors. T a is defined by 

Txx' ,yy' \zz' ~ ^^(^xX y X z ) ^ A xyz [m] 



-tH x 



A x , y , z ,[m] {\' x \' y \' z 



.1/2 



(12) 
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FIG. 2: Flow chart of the one-step renormalization of the 
wave function, (a) To use exp(— Hi^ +X r) to act on the ten- 
sor network state, (b) To evaluate the ^-matrix defined by 
Eq. ©. (c) Perform the singular value decomposition for S. 
(d) Truncate the basis space of X x and find A and B with Eqs. 
((8) and (|9]), respectively. 

T b is similarly defined. Thus we can apply the ten- 
sor renormalization group method proposed by Levin et 
al.[n| to evaluate and (^|6|*). 

To perform the tensor renormalization, we first take 
two T a and T b on the two ends of a bond and define the 
following D A x D A matrix 

M WM ,=Y,TZv k T b nkll . (13) 

n 

By taking the singular value decomposition, one can also 
express this matrix as 

M W , kk , = Uu>, n AnVkk',n, (14) 

n=l---£> 4 

where U and V are unitary matrices, A n is a positive 
defined diagonal matrix of dimension D A . Again we will 
truncate the basis space and keep only the basis states 
corresponding to the largest D 2 singular values of A. 
Then the M matrix can be approximately expressed as 

Mu> ik k' ~ Snll' S nkk>i 
n=l---D 2 

where 

S a nW = y/KUw, n , (16) 

Snkk' = V^nVkk^n (17) 

are the two vertex tensors defined in the new lattice 
shown in Fig. (j3^). 

After the above transformation, the lattice structure is 
changed (Fig. [3)3). Now we replace each smallest triangle 
by a single lattice point. This introduce a coarse-grained 
honeycomb lattice with two coarse-grained tensors T a 




(a) (b) 

FIG. 3: Steps of coarse-grain: (a) Form the M matrix by 
tracing out the common bond indices of tensors T a and T b 
defined on the two neighboring sites with Eq. JHJ), and then 
perform the singular value decomposition and find two new 
tensors S a and S b defined by Eqs. flS) and (fT7)) . (b) Trace 
out all common bond indices of S^-tensors (similarly for the 
S'k-tensors) on a triangle formed by the three closed vertices 
to form a coarse-grained tensor T a defined by Eq. JT5J). 




FIG. 4: Tensor renormalization transformation on the honey- 
comb lattice. 

and f b defined by 

J-xyz xik yji zkj V 10 / 

ijk 

rjib _ qb qb qb /-|Q\ 

xyz / _j u xik u yji u zkj ' K- 1 -^/ 

ijk 

This coarse grain transformation reduces the lattice by a 
factor of 3 at each iteration. Iterating this procedure, at 
the end the honeycomb lattice will eventually becomes 6 
(Fig. [4]). One can then trace out all bond indices to find 
the norm of the wavefunction. 

The above coarse grain tensor renormalization group 
transformation can be straightforwardly extended to 
evaluate (*|6|*). The difference is that T a and T b 
now may become site dependent and their definitions are 
changed. 

We have applied the above approach the spin-^ anti- 
ferromagnetic Heisenberg model (j2]). Both the ground 
state energy and the staggered magnetization M defined 

by 

M(h) = (20) 

are calculated. In Eq. (|2Q|) . E(h) is the ground state 
energy in a finite staggered magnetic field h. The lattice 
size is N = 6 x 3 10 . The finite size effect is negligible 
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FIG. 5: (color online) The staggered magnetization M(h) as 
a function of the staggered magnetic field, at different D. 



compared with the truncation error resulted from the 
coarse grain renormalization. 

Table [J shows the ground state energy and the stag- 
gered magnetization as a function of D for the Heisenberg 
model with h = 0. The zero field staggered magnetiza- 
tion is obtained by extrapolating M(h) obtained at finite 
h (Fig. ED to the limit h 0. With D = 8, we find that 
the ground state energy E = —0.5506 and the staggered 
magnetization M = 0.21 ± 0.01 in the zero field limit. 
They agree well with the results obtained by the other 
approaches (see Table HI]) . 

In conclusion, we have proposed a novel method to 
treat the tensor-network wave function of quantum lat- 
tice models in two dimensions. It allows us to treat a 
tensor-network state with D as high as 8 accurately and 
efficiently. The ground state energy and the staggered 
magnetization of the S=l/2 Heisenberg model on the 
honeycomb lattice obtained with this method are con- 
sistent with those obtained by other methods. 

We acknowledge the support of NSF-China and the 
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TABLE I: The ground state energy per site E and the stag- 
gered magnetization M in the zero field limit as a function of 
D. 



D 


E 


M 


3 


-0.5365 


0.249 


4 


-0.5456 


0.228 


5 


-0.5488 


0.220 


6 


-0.5513 


0.206 


7 


-0.5490 


0.216 


8 


-0.5506 


0.212 



TABLE II: Comparison of our results with those obtained by 
the other approaches for the ground state energy per site E 
and the staggered magnetization M of the Heisenberg model 
with h = 0. 
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Method 


E 


M 


Spin wave [12] 


-0.5489 


0.24 


Series expansion [13] 


-0.5443 


0.27 


Monte Carlo [14] 


-0.5450 


0.22 


Ours D=8 


-0.5506 


0.21 ± 0.01 



